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ABSTRACT 

We present an independent estimate of the interstellar extinction law for the 
Spitzer IRAC bands as well as a first attempt at extending the law to the 24/im 
MIPS band. The source data for these measurements are observations of five 
nearby star-forming regions: the Orion A cloud, NGC 2068/71, NGC 2024/23, 
Serpens and Ophiuchus. Color excess ratios Eh-kJ Ek^-[\] were measured for 
stars without infrared excess dust emission from circumstellar disks/envelopes. 
For four of these five regions, the extinction laws are similar at all wavelengths 
and differ systematically from a previous determination of the extinction law, 
which was dominated by the diffuse ISM, derived for the IRAC bands. This 
difference could be due to the difference in the dust properties of the dense 
molecular clouds observed here and those of the diffuse ISM. The extinction 
law at longer wavelengths toward the Ophiuchus region lies between that to the 
other four regions studied here and that for the ISM. In addition, we extended 
our extinction law determination to 24/xm for Serpens and NGC 2068/71 using 
Spitzer MIPS data. We compare these results against several ISO extinction law 
determinations, although in each case there are assumptions which make absolute 
comparison uncertain. However, our work confirms a relatively flatter extinction 
curve from 4 - 8/im than the previously assumed standard, as noted by all of 
these recent studies. The extinction law at 24/im is consistent with previous 
measurements and models, although there are relatively large uncertainties. 
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Introduction 



With the launch of the Spitzer Space Telescope, identification of embedded young stellar 
objects (YSOs) from their infrared excess emission at mid-infrared wavelengths (3 — Sfim 
using IRAC, the Infrared Array Camera, and at 24/im, using MIPS, the Multiband Imaging 
Photometer) has been facilitated by Spitzer 's unprecedented sensitivity. These wavelengths 
are sensitive to the emission from disk and envelope material surrounding YSOs, making it 
easier to differentiate young stellar objects from foreg roun d or background stars where the 
emission is primarily photospheric. lAllen et al.l (120041 ) and iMegeath et al.l (120041 ) develop a 
classification scheme t o differentiate young stellar objects based on their IRAC colors, and 
MuzeroUe et al.l (120041 ) present a scheme utilizing 24/im data in combination with IRAC data. 



Gutermuth J2005h has recently presented a robust classification scheme for YSOs which uti- 



lizes de-reddened mid-IR fiuxes of the cluster members in combination with ground-based 
ph otometry. These cla s sification schemes are briefly compared in the recent review article 
by I Allen et al.l (l2007al ). iLada et al.l (120061 ) also employ de- reddened IRAC observations to 
identify YSOs in IC348, using an alternate classiflcation scheme. Since these schemes in- 
volve de-reddening, it is critical to provide an accurate extinction law at IRAC and MIPS 
wavelengths that is appropriate for young clusters embedded in nascent molecular cloud 
material. 

Examining the extinction to background stars beyond these embedded young clusters 
probes the column density of dust in the molecular cloud. The availability of 2MASS has 
allowed for construction of large scale extincton maps in the near-infrared wh ich trace a 
substantial range in column densities (ILada et al.l 120071 : iLombardi et al.l l2006l ) . With an 
appropriate extinction law, this technique could be extended to longer wavelengths where 
the extinction is lower and higher column densities can be probed. 



Throughout much of the IRAC range of interest (0.7 - 5/im). iDraind (119891 ) fit data ex- 
isting at that time to a power-law Ax oc A~'^, with (3 = 1.75, and extrapolated the power law 
to meet the beginning of a silicat e contribution at ~7 /xm in his F ig. 1; other early papers 
have found rangi ng from 1.61 jRieke fc Lebofskvlll98 j) to 1 80 jMartin fc WhittetlE990 



Whittet et al.lll993l ). Dust models (jWeingartner &: Drainell200ll ) employing size distributions 
of various grain components with known optical properties have been constructed to match 
observed extinction laws. Measurements of the extinction law from ~ 0.9/im to ~ 2/im 
have been reported by many investigators using a variety of techniques (from both ground- 
based and space observations) along a number of lines of sig ht, and these were thought 



to be "universal", i.e. generally in agreement with each other (IDraind Il989l: ICardelli et al. 



2001 



1989 ) , although recently some authors have questioned this universality (iGomez fc Kenyon 



Racca et al.ll2002l ). Longward of 2/im the extinction is less certain (see the range in 
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j3 above), and the most robust estimates are generally obtained with instruments on space- 
borne telescopes. Pre-Spitzer studies have relied primarily oi i measurenaents of hydrogen 
recombination or hydrogen rotationa . 
Rosenthal. Bertoldi fc Drapatzl 200ol: 



and ro-vibrational lines (iLutz et al.l Il996l : iLuta Il999 



Bertoldi et al.lll999l ) to estimate the extinction law for 



A > 2.5/im. For example, Lutzl ( 19991 ) used hydrogen recombination line emission toward 
the galactic center to probe the extinction from 2.5/im < A < 20/im and found that the 
extinction law Ax/ Ay is flatter from 4 — 8/im than the power laws not ed above, a nd ap- 
proximately 4x smaller than the value of the standard curve at 7/im. iLutzl (119991 ) notes 
that the lack of a deep minimum near 7/im is i n conflict wi th predictions as s uming the 
standard mixture of silicate and grap hite grains (jDraindll989l ). iBertoldi et al.l (119991 ) and 
Rosenthal. Bertoldi fc Drapatzl J2000h observed H2 rotational and ro-vibrational lines toward 



OMC-1 and assumed a power law (A~^'^ extending from 2.4 to 6fim), and fit two silicate 
features to the extinction law beyond 8yum. The latter authors also determined extinctions 
from hydrogen recombination lines, and obtained smaller values of the relative extinction by 
this method than with H2 lines. In that case A\/Av was also relatively flatter from 3/im - 
7/im, an d from 1 . 5-2 t imes smaller than the comparison curve (in this case the IR extinction 
law from ICohenI (119931 ) with modifications to the 1 0/xm sil i cate f eature strength: the Cohen 
law is very similar otherwise to the one adopted by iDraind (1l989l )). A study using ISOGAL, 
D ENIS and 2MASS data of intermediate AGB and RGB stars in the inner galactic plane 
bv I Jiang et al. I J2006h determined values for A^/Ak, and Ai^/Ak, that are more comparable 
to those of Lutz (11999) tha n to earlier 'standard' interstellar extinction curves in the solar 
neighborhood ( Drainel 1989 ). 

Using observations of stars they identified as reddene d clump giant star s in the galactic 
plane obtained as part of the Spitzer GLIMPSE survey, llndebetouw et al.l (120051 ) directly 
measured the extinction law for the Spitzer IRAC bands at 3.6, 4.5, 5.8, and 8.0/tm. They 
observed the massive star-forming region RCW 49 at a distance exceeding 2 kpc, as well as 
two off-field regions very close to the galactic plane, and supplemented the Spitzer data with 
2MASS observations of the same fields. For the off-field regions, their method determines 
the diffuse interstellar medium extinction, assuming that the red clump giants are uniformly 
distributed along the line of sight. Reddened red clump giant stars were selected from their 
position in the J vs. J — Kg color-magnitude diagram, and a reddening band was identified 
which was then used to determine Aj/Ak, assuming a narrow distribution of intrinsic colors 
and luminosities of the clump giants. Ax/Ak^ was then determined from [A] — Kg vs. J — Kg 
color-color diagrams for these objects. The extinction curve for all three of these regions was 
nearly identical and also appeared to have a shallower slope for A > 4/im than the standard 
law. R. Indebetouw kindly provided the data for these three regions so that we could process 
them using the same methods that we used to process our data, hence making possible a 
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direct comparison of the extincti on laws toward the youn g clusters studied here, and that 
along the lines-of-sight probed by llndebetouw et al.l (120051 ) . 



Here we determine the extinction law from 3 — 8yum toward five young, nearby (< 1 kpc) 
cluster regions, including the Orion A cloud, NGC 2068/71, NGC 2024/23, Ophiuchus and 
Serpens. We further determine the extinction at 24/im for two of the regions, NGC 2068/71 
and Serpens. In all cases we measure extinctions in relatively broad wavebands as determined 
by fixed, cold filters. By examining those stars without disk/envelope infrared excess emission 
in the cloud direction, we evaluate the selective extinction at mid-infrared wavelengths: 
the extinction derives primarily from their host molecular clouds, although extinction from 
the diffuse interstellar medium is also a factor. We compare our determination against 
prior determ inations in th e literature including the "standard extinction law" such as that 
compiled by iDraind (119891 ). It is important when comparing these determinations to keep 
in mind that spectroscopic measurements may differ from our broad-band measurements at 
these wavelengths. 

Our separate extinction law study at the wavelengths of the IRAC bands is motivated 
by a discrepancy between slopes of the bands of reddened sources in the H — Kg vs. Kg — [A] 
color-color diagrams for our clusters, and slopes found using the GLIMPSE data. Reddened 
non-excess sources form a band on the H — Kg vs. Kg — [A] color-color diagrams extending 
from the intrinsic colors of the stars parallel to the reddening vector whose slope is the 
color-excess ratio Eh-KsI Eks-\\\ = {^h — Aks)/{^Ks ~ ^x)- These color excess ratios, 
derived in Section [2] and given in Table 2 for our data, and for the GLIMPSE data provided 
by Indebetouw, illustrate at a glance the d i screp ancy between selective extinction slopes 
toward our sources and theirs. iHuard et al.l (120071 ) also found a slope discrepancy between 
the selective extinction law toward i solated dense mol ecular cores observed in the Spitzer 
Co res to Disks (c2d) Legacy program (lEvans et al.ll2003l ) , and toward those regions examined 
by llndebetouw et al.l (120051 ). 



For these young clusters, and for the regions examined by llndebetouw et al.l (120051 ) . 
we utilize in our analysis th ose stars without a mid-infrared excess in the IRAC color-color 
diagrams (jAUen et al.l 120041 ) . namely class III cluster objects, and field stars. Background 
dwarfs and red giants, including the red clump giants, are the primary probes of the extinc- 
tion that originate in the clusters' nascent molecular clouds. Using complementary 2MASS 
ground-based H and Kg data as well as Spitzer data on these objects, we measure the color 
excess ratios Eh_kJEk^^[x\ from H — Kg vs. Kg — [A] color-color diagrams, and from these 
ultimately derive A\/Ak^ for an assumed value of Ah/ A^^- The same method allows exten- 
sion to 24/im for the Serpens and NGC 2068/71 clusters, although these results are more 
uncertain because of the comparative paucity of detected sources. 
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In Section [2] the observations and data reduction will be discussed, and Section [3] will de- 
scribe source selection and the derivation of the reddening law from these data. SectionlHwill 
report the results and discuss the differen ces between the extinctio n law derived here for our 
data and for the data used in the study by llndebetouw et al. J2005h. as well as the extinction 

law derived using different methods and extinction probes by Indebetouw et ali (2005), with 

further reference to other determinations (Lutz 19991 : Jiang et al. 20061 : Rosenthal. Bertoldi fc Drapatz 



2000l : iHuard et al.l 120071 ). We have not determined the total to selective extinction in this 



study, so shall concentrate on relative extinctions, assuming a value for Ah/Ak, in order 
to make comparisons with other work. Because a non-universal extinction law may lead to 
variations in this assumed value, we consider how this assumption affects our results. 



2. Observations and Analysis 



Observations of several young star-forming clusters were obtained with the IRAC camera 
( iFazio et al.ll2004l ) aboard Spitzer as part of several surveys concentrating on star-formation 
within 1 kpc of the Sun, including the "Young Cluster Survey" (PID 6), the "Orion Survey" 
(PID 43) from G. Fazio's guaranteed observation time, as well as the "Cores to Disks Legacy 
Program" (PID 174 and 177). The IRAC camera takes simultaneous images in f our channels 



whos e isophotal wavelengths are 3.550, 4.439, 5.731 and 7.872/im respectively (IFazio et al. 



20041 ). Some of the regions observed cover several square degrees so that multiple AORs 
with the IRAC camera were required. NGC 2024 and NGC 2023 were imaged as part of the 
same region, while NGC 2068 and NGC 2071 were also observed together. Observations of 
the Or ion A cloud inclu ded coverage of OMC 2/3, ONC and L1641. High dynamic range 
mode (IFazio et al.l |2004| ) was used, resulting in a frame with a 0.4 second exposure time 
and a frame with a 10.4 second exposure time at each position. The maps for each of 
NGC 2024/23, NGC 2068/71 and Orion A were repeated twice. The maps of Serpens and 
Ophiuchus were repeated four times, with small offsets. The additional dithers provide both 
additional sensitivity to faint sources and greater ability to remove cosmic rays from the final 
mosaics. We have not corrected for the position dependent gain variation for point source^, 
which would result in an additional 5% uncertainty. 

For the Orion fields individual frames were mosaiced using a custom IDL program and 
from these mosaics, source selection and point sou rce photometry were completed using 
Gutermuth's Phot Vis vl.08 (iGutermuth et al.ll2004l ). The photometry used an aperture of 
2.4" and a sky annulus with inner and outer radii of 2.4" and 7.2". Details on the data 



http : //ssc . spitzer . caltech. edu/irac/locationcolor/ 
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reduct ion a nd source s e lection i n Serpens and Oph iuchus are described in IWinston et al. 
fl2007h and kUen et~aD J2007bh . iReach et all J2005h discusses the calibration of the IRAC 
instrument. The calibration of the Spitzer telescope was refined during the course of this 
work, and the resulting zero points are listed in Table 1. However, changes in the zero 
point do not change the results in this paper and are taken into account in our fitting 
technique. From the 2MASS point source database we obtained corresponding J, H and Kg 
band magnitudes for each star. For the NGC 2068/71, NGC 2024/23 and Serpens regions, 
we required magnitude uncertainty a < 0.1 in all seven bands for sources used to calculate 
the extinction law. For the Orion A cloud and Ophiuchus a detection in J band was not 
required, thereby increasing the number of highly reddened sources and providing a better 
determination of the extinction law. For the remaining 6 bands, a magnitude uncertainty 
of (J < 0.1 was imposed. If we do not require a J band detection for NGC 2024/23 , NGC 
2068/71 or Serpens, our derived values for the slope of the extinction vector are within la 
of the values derived with the seven-band data. 

Requiring sources to have this magnitude uncertainty reduced the spread of the sources 
on each color-color diagram, and minimized the uncertainty in the fit (see Section [3]). Those 
sources for which only upper limits to 2MASS magnitudes were available were not included 
in our analysis. Most of the sources had much smaller magnitude uncertainties. The median 
magnitude uncertainties for NGC 2024/23 sources are 0.03, 0.02, 0.02, 0.005, 0.005, 0.02, 0.05 
for J, H, Kg, and channels 1 through 4 of IRAC respectively. Similar median magnitude 
uncertainties were obtained for the other regions studied here. Loosening the requirement 
on the maximum uncertainty will include more sources with lower fluxes. It is possible 
that these sources have higher extinction and could provide a more accurate measure of 
the reddening law. However, increasing the uncertainty limit leads to the inclusion of more 
outliers that skew the slope and are difficult to remove with a simple color cut. A small 
number of fits did not have any outliers and for these fits including the sources with higher 
uncertainty does not change our results by more than the measured uncertainty in the color 
excess ratios. 

MIPS 24 /im images were available for several of the regions studied. Not all regions 
were suitable - there were insufficient non-excess sources in some regions to reliably derive 
an extinction law. We report on NGC 2068/ 71 and Serpens o bservations with the Multi- 
band Imaging Photometer for Spitzer (MIPS, iRieke et al.l 120041 ). We summarize the MIPS 
observations, with more detailed descripti ons of the data reduction and photometri c mea- 
surements for Serpens and NGC 2068/71 in I Winston et al.l (120071 ) and lMuzeroUe et al.l (120071 ) 
respectively. The measurements were obtained in scan mode with half-array cross-section 
offsets. All three MIPS bands are observed simultaneously in this mode; the total effective 
exposure time per point is about 80 seconds at 24 microns. Photometry was determined 
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using a PSF fitting program in IDL. 



3. Background Source Selection and Determination of Extinction Law 



Background dwarfs and giants (including red clump giants and RGB stars) as well 
as embedded pre-main sequence stars can be used to measure the extinction due to the 
molecular cloud dust associated with the star-forming region as well as diffuse line-of-sight 
extinction. As mentioned above, the sources used to measure the extinction were chosen 
from those without infrared excesses from disk or envelope emission in the [3.6] — [4.5] vs. 
[5.8] — [8.0] color-color diagram toward each cluster region. Eliminating infrared excess 
sources from consideration removes any ambiguities between stars with red colors caused 
by extinction vs. stars with circumstellar dust or envelope emission, including both pre- 
main sequence stars and AGB stars. A non-excess source was defined as having colors 
[3.6] — [4.5] < 0.6 and [5.8] — [8.0] < 0.2. This region of consideration was chosen so that it 
di d not overlap with the expected position of pre-main sequence stars with disks according 
to lAUen et al.l (120041 ) and included the expected spread in colors due to interstellar dust 
reddening, which is larger at the shorter wa velengths. This definition of non-excess sources 
is consistent w i th the the oretical models of IWhitnev et al.l (120031 ) . and the observations of 
Megeath et al.l (120041 ) and iGutermuth et al.l (120041 ). An upper limit on the [3.6] — [4.5] color 
is used to eliminate unresolved outfiow sources misidentified as stars. Outfiows have emission 
features in the 4.5 /im band, leading to [3.6] — [4.5] colors which could be confused with those 
from reddened non-excess stars. 

We extend our selection criteria to attempt to further remove foreground stars, and 
stars whose colors place them outside the reddening band. We re move foreground giants 



by ex cluding any star brighter than a K2III giant with Mj = —0.9 (iLopez-Corredoira et al. 



20021 ) at the distance of each cluster. We also exclude stars with H — Kg < 0.2, which are 
mostly foreground dwarfs (see discussion below). Stars with Kg — [3.6] < do not fit into 
the reddening band and are removed as well as stars that appear to have an excess in the 
H~KsVS i^s- [8.0] diagram (e.g. H-Ks< 1.04fK, - [8.0]) -0.2 for Serpens). We have not 
attempted to use only clump red giants, as was done by llndebetouw et al.l (120051 ) . because of 
the difficulty in distinguishing these stars from highly reddened dwarfs and pre-main sequence 
stars without disks. Our reddest stars will be dominated by red giants and embedded pre- 
main sequence stars because of the large intrinsic luminosity needed to detect these stars at 
such high extinction although at lower extinction we may have a mix of reddened dwarfs as 
well as giants and pre-main sequence stars. For Orion A, a number of galaxies contaminated 
the field and their removal from consideration is discussed in section 13.11 
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Lopez-Corredoira et al.l (|2002| ) point out that stars lying to bluer colors than the red 
clump giant distribution in a J or vs. J — K color-magnitude diagram for a given line of 
sight are primarily dwarfs, while stars lying to redder colors are predominantly M giants and 
AGB stars. The red clump giant distribution is almost vertical around a median J — K = 
0.75 ± 0.2 for off-cloud lines of sight at galactic latitudes with absolute value > 6°. For fields 
probing greater extinction from either molecular clouds or along lines of sight closer to the 
galactic plane, the red clump giant distribution shifts to redder J — K, fainter magnitudes 
and the distribution is broader. All of our clusters have galactic latitudes exceeding 6°, with 
the exception of Serpens, at 6=5.4°. To illustrate the above points for one of our clusters, we 
show a J vs. J — Kg diagram for NGC 2068/71 (a region of high extinction) as well as an off- 
cloud region overplotted (Fig. [1]). The off-cloud field is centered at 1=210°, b=-8°, a region 
that was chosen because it did not show any significant CO emission indicative of a dense 



molecular cloud (IWilson et al.ll2005l ). In that diagram, the sources coincident with the cluster 
utilized in measuring the color-excess ratios are plotted with blue symbols, and red symbols 
illustrate cluster sources rejected by the H — Kg cutoff. Off-cluster sources are plotted with 
black symbols. In the black double-peaked distribution, dwarfs at various distances along 
the line of sight occupy the bluer peak, and the red clump sources the redder, with median 
color (J — K)o = 0.75, as noted above. The effect of the molecular cloud's reddening can 
be seen by noting that most red clump extinction probes to NGC 2068/71 form a clump 
centered at J — Ks=1.05 at J=14.5. The off-cloud giants have nearly constant J — Kg color, 
indicating little extinction from the diffuse ISM. For this reason we expect the contribution 
of the diffuse ISM to our measured extinction law for the Orion fields to be small, although 
this may not be equally valid for Serpens and Ophiuchus. Many sources redder than the red 
clump giants are sampled, and in particular, for the cluster region (red and blue points) in 
Figure [1] there is no distinct red clump giant locus. In the cluster region, the highly variable 
extinction makes it impossible to isolate red clump giants, reddened dwarfs and embedded 
pre-main sequence stars. We are likely dominated by red clump giants, and perhaps pre- 
main sequence stars, in the densest parts of the cluster, because they are bright enough to 
be observed at large reddening. However, we do not attempt to use only red clump giants. 

We now move to discussion of reddening law determination. A straight line was fit to 
background sources in the 'reddening' band on an if — Kg vs. Kg — [A] color-color diagram. 
After all of the restrictions described above were imposed, there were 415, 670, 1830 and 
1225 background sources for NGC 2023/24, NGC 2068/71, Orion A and Serpens respectively. 
Ophiuchus had 120 background sources and the implications of this small sample size are 
discussed further in section HI The slope of the hne fit to the 'reddening' band gives the 
color-excess ratio Eh-kJ Ek^-[\] which is then used to determine the rela tive extinction law 



A\/Aks, where we assume Ah/Ak^ = 1.55 ± 0.08 (the value derived by llndebetouw et al. 
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( 120051 ) )■ The effect of a different choice of Ah/Ak^ is discussed in section I4.1.3I 



The hne fit to the reddening band data in each of the H — Kg vs. Kg — [A] color-color 
diagrams minimized the chi-squared, weighted by the uncertainty in each color. 

Xia,b) = ^ (1) 

where a and h are the y-intercept and slope, respectively, of the best fit line and a^i and 
ayi are the uncertainties in the x and y-values respectively. The algorithm amoeba in IDL 
was used to minimize the chi-squared. Determining the best fit-line to the reddening band 
data using equation [T] not only accounts for the magnitude uncertainties, but also provides 
a method for estimating the uncertainty in the slope, which leads directly to A\/Ak^i the 
quantity of interest. 



Press et al.l (119921 ) presented a detailed analysis of the adopted fitting technique as well 
as the characterization of its uncertainty, which we summarize here for the convenience of the 
reader. At the best-fit points the is a minimum, but if the values of a and b are perturbed 
away from the best-fit values increases. A constant value of Ax^ = — Xmin gives a 
confidence region within parameter space, which depending on the size of Ax^, could enclose 
a 68%, 90%, etc. confidence interval. The uncertainty in a and b is then the projection of 
this confidence region onto the a and b axis. Finding a, b where Ax^ has a specific value 
is determined analytically with a Taylor expansion when is linear in 6, but must be 
determined numerically here because x^ is non-linear in b. This determination is performed 
by varying 6, minimizing a at each step, until Ax^ equals the desired value. There are two 
values of 6, one greater and one less than the best-fit value, for which Ax^ equals the desired 
value. The distances of these two new values of b from the best fit slope are combined to 
determine the uncertainty in the slope. 

For the majority of our analysis Ax^ was chosen to be 1, which corresponds to a 



68% confidence level (jPress et al.lll992l ). Increasing the confidence level to 99.99% requires 



Ax = 15.1, which results in an increased uncertainty in the slope. Both a 68% and 99.99% 



confidence interval were determined and are included in Table |5J The uncertainties are ~ 4 
times bigger for the larger confidence interval than for the smaller confidence interval, but 
lead to at the most a 6% uncertainty in the fit. 

The process that was used to derive the extinction law for the IRAC bands can be applied 
to the 24/im channel of MIPS. There were only sufficient numbers of non-excess sources with 
corresponding detections in 2MASS, IRAC and MIPS in the Serpens and NGC 2068/71 
regions. The non-excess sources were defined as having [3.6] — [4.5] < 0.6, [5.8] — [8.0] < 0.2 
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and [8.0] — [24] < 0.5, with the last color cut used to remove sources that had an excess only 
at 24yum. For both NGC 2068/71 and Serpens fewer than 50 sources were used for the final 
fit, making sure to exclude outliers and sources with if — < 0.2, so we are working with 
quite small numbers of sources. For this reason, we view these results with extreme caution. 

3.1. Orion A 

While the process for deriving the extinction law from Serpens, Ophiuchus, NGC 2068/71 
and NGC 2024/23 data was straightforward, the number and variety of sources toward Orion 
A required more careful attention. The Orion A cloud covers 4.4 square degrees, a much 
larger area than the other four regions observed, increasing the number of sources used when 
determining the reddening law but also increasing the number of outliers. These outliers are 
most obvious on the H — Kg vs. Kg — [8.0] color-color diagram on which there is a 'branch' 
with H — Kg ^ 0.25 which could skew the fit. 

This 'branch' is not as dominant in any of the other H — Kg vs. Kg — [A] color-color- 
diagrams and has nearly constant H — Kg c:^ 0.25 color as well as J — Kg ~ 0.9. The 
objects in this branch are too blue in H — Kg color to be infrared excess sources but they 
are too red in Kg — [8.0] color to be main-sequence stars, red giants, or Class III sources. 
The most hkely explanation is that these sources are galaxies. Their blue colors in the near 
infrared likely derive from main-sequence stars dominating the light at short wavelengths. 
Polycyclic Aromatic Hydrocarbon (PAH) emission from the dust within the galaxies is the 
likely cause of the red Kg — [8.0] color, since strong PAH feature emission falls within the 
[8.0] band of IRAC. Interpreting these sources as galaxies is consistent with their fairly 
uniform spread in position throughout the Orion A cloud. They arc also not present in 
the areas with the highest concentration of young stellar objects because these areas have 
a large density of nebulosity which decreases the sensitivity for detection of background 
objects, and a high column density of dust obscuring faint background sources. It is also 
possible that foreground sources contaminated by extended PAH emission are contaminants: 
however, this does not explain the lack of outhers in the regions of largest PAH emission 
near the young stars. Regardless of the physical nature of these outliers, we remove them 
from our fitting procedure. A first attempt at removing these outliers was performed by 
changing the selection criteria for non-excess sources. This did not effectively remove the 
outliers without also removing a signifcant number of sources in the reddening band which 
are needed to accurately measure the extinction law. Removing these outliers, without 
removing a substantial number of sources in the reddening band, was completed by first 
fitting a fine in the H — Kg vs. Kg — [8.0] color-color diagram to all of the non-excess 
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sources, including the outliers, and then using the slope of this fit to define a boundary line 
{H — Kg = l.ll{Ks — [8.0]) — 0.2), below which were assumed outliers and above which were 
the sources used for a more accurate measure of the extinction. The y-intercept of this line 
was chosen so that a majority of the outliers was removed without removing a significant 
portion of the reddening band. The initial fit in the H — Kg vs. Kg — [8.0] color-color diagram 
had a slope of 1.110 and the final fit had a slope of 1.116 with x^/^ decreasing from 1.99 
for the initial fit to 1.48 for the final fit. While this is a small correction, the final fit is a 
more accurate measurement of the extinction law at 8/im. 



4. Results 
4.1. The Extinction Law from 3-8/im 

4.1.1. IRAC-Derived Color Excess Ratios 

The slopes of the linear fits to the H — Kg vs. Kg — [A] color-color diagrams for our five 
clusters are listed in Table [2] and plotted in Figures El [3], HI [5], and [61 It is apparent that 
our derived slopes (color excess ratios) are reasonably consistent from region to region, but 
in a few instances, are more than 1-a different from each other at the longer wavelengths. In 
particular, the selective color excess ratio toward the Ophiuchus cloud is lower than toward 
the other four clouds at 8.0/im. Since Ophiuchus has the smallest number of non-excess 
sources (120 vs. >400 for the other four regions) we assess the possibility that the difference 
is due to the small number of non-excess sources. We randomly select 120 non-excess sources 
from the Orion A field and measure the color-excess ratios in the four IRAC bands for the 
subsample, repeating this process for 1000 subsamples. We can estimate how likely it is to 
measure the Ophiuchus color-excess ratios based on the distribution of color-excess ratios in 
each of the four bands for the subsamples. For the 5.8 and 8/im bands the average color- 
excess ratio for the subsamples is 1.115 and 1.122, with a standard deviation of 0.06 for 
both bands. The Ophiuchus color-excess ratios are approximately 2a from these mean color- 
excess ratios at both 5.8 and 8/im, only a marginally significant deviation. The fraction of 
subsamples with color-excess ratios at or below the Ophiuchus values are 3% and 1% for the 
5.8 and 8/im bands respectively. Thus sample size is unlikely to be the cause of the lower 
selective color excess ratio at 8.0/im for Ophiuchus. 

The color excess ratios for the / = 284°, 6 = 0.25° off-cloud line-of-sight GLIMPSE data 
were derived using the same source selection criteria and methods as used for the five star- 
forming regions discussed here. In contrast to the general agreement among the clusters, 
the slopes from the GLIMPSE data listed in Table [2l are generally different from the cluster 
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slopes. For completeness, we also quote the color excess ratios from llndebetouw et al.l (120051 ) 
for the / = 284° off-cloud line of sight. Since they directly measured E[x]-Ks/ Ej-Ks^ we form 
Eh-KsI Eks-[\] from those ratios in Tabled Those values differ at 3.6 and 4.5 /im from 
values determined by our analysis. 



The question remains whether the larger sample size employed by llndebetouw et al. 



( I2OO5I ) in determining the extinction law is the cause of the discrepancy, as was suspected 
for Ophiuchus. Their method uses red clump giant stars selected by their J — Kg colors and 
only requires a detection in the three bands {J, Kg, [A]) that are used when measuring the 
color-excess ratios. This means that more sources are available for determining the color- 
excess ratios. To av oid this, and any o t her p ossible biases due to differing methodologies, we 
have reanalyzed the llndebetouw et al. J2005h data using the same methods that were applied 
to our five clusters as noted above. The I = 284°, b = 0.25° off-cloud line-of-sight was studied 
in detail because it had the largest number of background sources and its extinction law was 
representative of the other two lines of sight. There were 1590 background sources that 
met our criteria, comparab l e to th e number of background sources used for four of our five 
clusters, llndebetouw et al.l (120051 ) find 3560 red clump giants along this same line of sight 
when requiring a detection in J and Kg band. This includes stars that have an infrared excess 
due to circums tellar dust and sta r s tha t are not detected by IRAC. The actual number of 
sources used by llndebetouw et al.l (120051 ) when measuring the color-excess ratios was smaller 
(< 3000) after requiring an IRAC detection and removing stars with an infrared excess. We 
examine the effect of differing sample sizes by considering the Orion A field and varying the 
limit on the magnitude uncertainty. When limiting the magnitude uncertainty to 0.1, there 
were 1825 background sources, while changing the limit from 0.05 to 0.25 results in sample 
sizes ranging from 1265 to 2360. The color-excess ratios for the extreme cases differ from 
the case where the magnitude uncertainty is less than 0.1, and from each other, by less than 
the measured uncertainty in the color-excess ratios. As a result we do not expect differing 
sample sizes to provide a source of systematic error when comparing our results to those of 
Indebetouw et al.1 (l2005l ). 



Figure [7] gives the color excess ratio Eh^Ks/ Eks-[X\ versus 1/A for the five regions 
observed here as well as for the GLIMPSE / = 284° off-cloud line-of-sight data analysed 
by us, showing the systematic separation between the color excess ratios derived for the 
GLIMPSE line of sight, and for the fiv e young clusters. Our va lues are consistently 1-a, or 
more, higher than those rederived from llndebetouw et al.l (120051 ) data with the exception of 
the 4.5/im band. 

The IRAC non-excess stars used when measuring the color-excess ratios can include class 
III pre-main sequence stars, main sequence stars and red giants (clump giants, and RGB 
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stars), all of which have different intrinsic colors. Our constraint to exclude unreddened stars 
and bright probes should minimize contribution from foreground objects. The proximity 
of Serpen s and Ophiuchus reduce s the number of foreground contaminants which must be 
removed. Ilndebetouw et al.l (120051 ) restrict themselves to red clump giant stars and find they 
are well-separated from main sequence stars on a J vs. J — Kg color-magnitude diagram for 
the regions they studied. 

We examine the dependence on extinction probe by calculating residuals which are de- 
fined as th e difference betwee n {H— K^)nhf; for a given i^^ — [A] and the predicted {H—Ks)ind05 
using the Ilndebetouw et al.l (120051 ) determination of Eh-kJ Ek^-\- Figure [H where the 
residuals are plotted versus H — Kg color, shows the deviation of our extinction law deter- 
mination for NGC 2068/71, using non-excess stars, from that of Ilndebetouw et al.l (120051) 
using red clump giants. As noted in § 3, the parameterization used by Ilndebetouw et al. 



( I2OO5I ) to pick out red clump giants will not apply to our cluster lines of sight. As can be 
seen in Figure [1] our extinction probes in NGC 2068/71 exhibit varying levels of extinction. 
If using only red clump giant probes produced a systematically different extinction law than 
we measured, we would expect a locus of sources to follow the Ilndebetouw et al.l (120051 ) law 
at all four bands. No such locus is observed in Figure [HI suggesting that extinction probes 
are not the main cause of the systematic difference. So urces with H — K,^ > 0.6 , which are 
mostly red giants, show a substantial deviation from the Ilndebetouw et al.l (120051 ) extinction 
law. Since both methods referred measurements to the i^'^-band rather than to the i^-band, 
a correction for different photometric systems is not necessary. 

Rederived GLIMPSE color excess ratios are systematically separated from the color 
excess ratios for the five young clusters observed here. While there are differences in the 
results derived from the GLIMPSE data using the different methodologies, especially at 3.6 
and 4.5/im, there is stil l a sig nificant separation between our color excess ratios and those 
from Ilndebetouw et al.l (120051 ). This suggests that the separation is not primarily caused 
by the differing methods and/or extinction probes. The systematic separation may instead 
be due to differences in the extinction from dust in the molecular clouds and the diffuse 
interstellar medium. This possibility will be further discussed in § 14.1.41 



4.1.2. IRAC-derived Extinction Law 



We derive the extinction law A\/Ak^ for the five clusters studied here as well as for the 
selected GLIMPSE data, using the color excess ratios determined above. Ax/Ak,, is related 
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to the color excess ratios, and the assumed value of Ah/Akq by: 

Ah 1 
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Since we have used precisely the same methods to establish the extinction law for both data 
sets, we have removed uncertainties in the comparison that are method-dependent. In both 
cases we assume Ah/Aks = 1-55: values for Ax/Ak^ are listed in TableO I n addition, we also 
give th e extinction law for the / = 284° off-cloud line of sight as derived by llndebetouw et al. 

Different choices for Ah/Ak^ will affect our derived values of Ax/Ak^- For NGC 
2024/23, for example, if Ah/Ak^ = 1-53 (/3 = 1.61) was chosen then Ax/Ak, would be 
0.64,0.55,0.51,0 52 for the four IRAC bands. If instead Ah/Ak, of 1.73 (/?=1.99) from 
Nishiyama et al.l (120061 ) was chosen then Ax/A^^ would be 0.51, 0.38, 0.32, 0.34 for the four 
IRAC bands. The difference in these two determinations of the relative extinction law is 
larger than the statistical uncertainties in Ax/Ak^ derived from the uncertainties in the se- 
lective color excess ratios, which are < 0.01. The uncertainties reported in Table [3] include 
the uncertainties in the measured color excess ratios but do not include the uncertainty in 
Ah/Ak,. 

The values of Ax/Ak^ for each of the five regions observed here, subject to the assump- 
tion of Ah/Ak, = 1.55, are plotted in Figure P The rederived extinction law for the I = 284° 



off-cloud line of sight frona llndebetouw et al.l (120051 ) is also plotted in Figure M Extinction 



law measurements by iLutzl (119991 ) , included as squares in Figure are more consistent with 
Indebetouw et all J2005h . a lthough sorn e of th eir measurments overlap with our results. The 



measurement at 7/im from iJiang et al.l (120061 ) is marked with a circle in Figure [9], and the 
spread in th eir results for ma n y diff erent sight lines encompass both our results and the 
results from llndebetouw et al.l (120051 ). Our extinction law deviates from a power law, as is 
also seen in these previous studies. 



4.1.3. Comparison with ISO-Derived Extinction Laws from 3 — 8jj,m 



Luta (119991 ) compared the observed to expected fluxes for 15 hydrogen recombination 



lines toward the galactic center from ISO data and from this determined the extinction ratios 
at the wavelengths of the emission lines. Theoretical estimates of the expected fluxes assumed 
Te = 7000K, ne=3000 cm~'^, case B recoi nbination (which they verif y from their observations) 
and refer directly to the calculations by IStorey Sz Hummer Jl995h . The recombination line 
results also depend on the assumed helium abundance, [HeII]/[HeIII], and the assumption 
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that line emissions emanate from the same volume. A crucial part of their analysis is 
the scaling of their relative extinction law: they deduce ^4^/^2.626 (Br/3 at 2.626/im is the 
shortest-wave emission line they observe) and use a A~^'^^ scaling law to refer to i^'-band. 
We will use a A~^-^^ law which corresponds to our assumed Ah/A^s = 1-55, allowing for a 
more direct comparison with our values of Ax/Ak^ (see discussion below). Numerical data 
on Ax/Av were provided courtesy of D. Lutz: here we deduce from these data Ax/ A2S26, and 
using a factor of A/^^/A2, 626=1-38 (assuming the power law index /5=1.65), deduce Ax/A^^ 
for their data. Since we use 2MASS JHKg data, we have taken the wavelength for f^s-band 
to be 2.16/im. Power law indices of 1.61 and 1.8 lead to at most a 2% error in the extinction 
ratio. T his is consid erably more certain than variations in Ax/Ay, which can be as much 
as 30% (lGlasslll999l ). so we choose not to compute Ax/Ay for our data. 



Jiang et al.l (120061 ) combined ISO observations at 7 and 15/im of the inner galactic plane 



with 2MASS and DENIS near-infrared observation to determine the mid-infrared extinction 
law and its variation along different sightlines. They assumed a constant value of (J — Ks)o, 
{Kg — [7])o and {Kg — [15])o for early AGB and RGB-tip stars and measured the color e xcess 
of these stars along various sightlines. Assuming Aj/Ak from lRieke fc Lebofskyi (1l985l ) they 
find At/Ak = 0.47±0.04 and Ai^/Ak = 0.40±0.07 (marked as circles in RgnieM)- Although 



Rieke fc Lebofskvl (I19851) d i d not measure the extinction law with the 2MASS photometric 



system, Indebetouw et al. ( 20051 ) found Aj/ Ak„ in the 2MASS photo metric system to be 
similar to the val ue of Aj/Ak measured by iRieke fc Lebofskyi ( 119851 ). The results from 



iang et al.l (120061 ) do not change significantly when using A,j/Ak„ from llndebetouw et al. 



(I2OO5I ) as opposed to Aj/A^ from iRieke fc Lebofskyi (119851 ). At 7yum they find the extinc- 



tion law measurements have a Gaussian distribu tion whose one sigrna siz e enc ompa s ses th e 
measurements presented here, as well as those of llndebetouw et al.l (120051 ) and iLutzl (119991 ). 
There is more uncertainty at 15/im but there also appears to be evidence for a Gaussian 
distribution of extinction values, Ai^/A^- 



We chose Ah/Ak, from llndebetouw et al.l (120051 ) to provide a direct comparison of all 
the mid-infrared selective exti r ictioii law measurements discussed here. The extinction law 
measurements bv iJiang et al.l (120061) use a similar normalization to our assumed Ah/Aks, 



and the measurements of iLutzl (119991 ) were normalized using a power law form that is con- 
sistent with Ah/Ak, = 1.55. Different assumptions of the near- infrared extinction law will 
affect our derived Ax/Ak^ as discussed earlier, but we choose a single consistent normaliza- 
tion to better compare the various determinations of the mid-infrared extinction law. 
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4- 1.4- Molecular Cloud vs. Diffuse Interstellar Extinction. 



A plausible explanat i on for the separation between our extinction law and the extinction 
law of llndebetouw et al.l (120051 ) is a variation in the extinction law between that through 
the diffuse ISM and through molecular clouds. Variation in dust optical propert ies between 
molecular cloud dust and that in the ISM has been invoked by lHuard et al.l (120071 ). who, from 
observations of cloud cores obtained as part of the c2d Spitzer Legacy p rogram, have also 



found extinction law slopes (color excess ratios) which deviate from those of llndebetouw et al. 



The cold ISM includes diffuse interstellar clouds {Ay < 0.1), poorly s hielded regions 



of mo lecular clouds {Ay < 3), and dense molecular cloud cores {Ay > 3). IWhittet et al. 
( 120011 ) and references therein, have explored grain growth in dense clouds, and have re- 
ported differing extinction laws for the diffuse ISM and dense clouds. For reference, Ay=3 
is s hown as a vertical dashe d line in Figure [HJ most of the deviation of our results from 
the llndebetouw et al.l (120051 ) extinction law occurs at higher extinctions. Changes in the 
extinction law originating in these regions are likely due to grain growth, either by coagu- 
lation of small grains or the growth of icy mantles, both of which will affect the continuum 
level of extinct i on. S mall grains are thought to coagulate to form larger grains (e.g. p Oph, 
Whittet et al.l (l200ll )): grain emissivity is dependent on the optical constituent properties, 
and for the same constituents is proportional to a power(powers) of a/ A where a is the grain 
size. Grain growth will not only remove small grains as potential absorbers, but also as 
emitters. The lack of e mission from small grains has been observed along lines of sight with 
high column densities (IStepnik et al.ll2003l : Idel Burgo et al.ll2003l : ICambresy et al.ll2005l ) sug- 
gesting grain growth le ading to porous grains with enhanced emissivity within these clouds. 
Pendleton et al.l (120061 ) have found a decrease in the depth of the 9.7 /im silicate feature 
with increasing Ay, indicating silicate grain growth. Ay = 3 is the approximate extinction 
threshold for the detection of ice m antles, and in the Taurus Dark Cloud, icy mantle growth 
is imphcated (I Whittet et al.l 1200 ll ). Dust in molecular cloud cores is protected both from the 
UV interstellar radiation field and from cosmic ray processin g; those dust grain s may there- 
fore have larger mantles consisting of water and organic ices (jGibb et al.ll2004l ). Absorption 
features between 3 and S/xm due to ices are often seen along lines of sight intersecting dense 
mole cular clouds (iKnez et al.l l2005l : iGibb et al.l |200J) while they are not seen in the diffuse 
ISM ( IWhittet et al.lll997l ). On the other hand, grain processing and fragmentation because 
of radiative shock propagation, and turbulence, may also be involved in grain evolution in 
both molecular clouds and in the diffuse ISM. 



Our extinction law for the young clusters s tudied here derives main ly from dust asso- 
ciated with their molecular clouds, while that of llndebetouw et al.l (120051 ) derives primarily 
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from diffuse interstellar dust, especially for the / = 284° off-cloud line of sight. There is 
a greater contribution from diffuse interstellar dust for the closest clusters in our sample, 
Ophiuchus and Serpens, partially because they are in the inner galaxy, and partially be- 
cause the extinction probes (not excluded by the JHK brightness cutoff) lie well behind 
the cluster, so that their light passes through more of the ISM before reaching the cluster. 
The extinction curve toward Ophiuchus has been t ermed anomalous: measurements of the 



UV/optical extinction curve and linear polarization (IGreen et al.lll992l : iKim fc Martirull996l ) 
toward Ophiuchus indicate a depletion of small grains and the presence of large grains within 
the molecular cloud. Studies of other dark clouds have produced similarly anomalous ex- 
tinction as compared with the ISM: we are unaware of any comprehensive study co mparing 



optical/UV extinction in the dark clouds studied here. iTeixeira fc EmersonI (119991 ) observe 
H2O and CO ices in the Taurus and Ophiuchus dark clouds, and concur with previous stud- 
ies that different physical conditions obtain in those two clouds. Although column densities 
of the H2O ices and the visual extinction are correlated in each case, the regression line is 
different. There are several Ophiuchus sources which follow the Taurus regression line rather 
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/? oc z/^) in the far-infrared and find variations between Orion, Ophiuchus and Serpens. 
This indicates that the grain properties differ between these three clouds, but the anticor- 
relation between spectral index and temperature within a cloud suggests that there can be 
significant variations in grain properties along different lines of sight through a single cloud. 
Consequently, we cannot determine here the extent to which the molecular cloud mid-IR 
extinction toward Ophiuchus is anomalous compared with other clouds in our sample. 

In contrast to our clusters, all within 500 pc of the Sun, the lines of sight for the ISO 
measurement may include a large r cont r ibutio n from a mixture of extinction from the diffuse 
ISM and molecular clouds. The iLutzl (119991 ) extinction law to the galactic center ionized 
gas is dominated by dust in the diffuse ISM along the line of sight, although sorn e of the 
extinction is due to molecular cloud material (iLutz et al.lll996l : IWhittet et al.l 119971 ). Broad 
ice absorption features at 3 and G/xm are seen towa rd Sgr A* an d are a ttributed to foreground 
molecular cloud material (jChiar et al.ll2000l ). The iJiang et al.l (120061 ) extinction law is for a 
mix of diffuse ISM and dense molecular cloud dust. Given the uncertainties in the derived 



laws of I Jiang et al.l (120061 ) and iLuta (119991 ) . we cann ot distinguish whether an y agreement, 
or lack of agreement with either our law or that of llndebetouw et al.l (120051 ) depends on 
whether the material causing extinction is part of the diffuse ISM or a dense cloud. 

The differing dust properties between diffuse ISM and dense molecular clouds could 
lead to differences in the extinction law. We directly measure Eh-KsI Eks-[\]^ which relates 
A\\t^/Ak, and the assumed value of Ah/Aks- A variation in the near-infrared extinction law. 
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the mid-infrared extinction law, or both, could lead to the observed systematic separation 
in the selective co l or ex cess ratios between our observations through molecular clouds and 
Indebetouw et al.l (120051 ) 's observations through the diffuse ISM. If a larger value of Ah/Ak^ 
was more appropriate for the extinction law toward molecular clouds, th en A^/Ak, would 



be sm aller and the mid-infrared portion of the extinction law deduced by llndebetouw et al. 



( I2OO5I ) and us could be made to agree. Values of Ah/A^^ between 1.6 and 1.65, depending 
on the region, would bring our mid-infrared extinct ion law into agreement with the average 
extinc tion law reported bv llndebetouw et al.l (l2005l ). assuming Ah/Ak^ = 1.55 still apphed 
to the Indebetouw et al. J2005h extinction law. It may also be the case that the near-infrared 
extinction law does not appreciably vary between the diffuse ISM and molecular clouds and 
there is instead a variation in the mid-infrared extinction law. We now examine the near- 
and mid-infrared data separately in order to attempt to break this degeneracy. 

The possibility of var iations in Ah/Ak, i s sugg ested by recent observations. Although 
Cardelli et al.l (119891 ) and iMartin fc Whittetl (jl990l ) found the near-infrared extinction law 
to be universal, more recent studies have produced mixed results. One of the difficulties 
in evaluating the variation in the near-infrared extinction law along different lines of sight 
is that different photome tric systems (e.g. K vs. K,^) can have a substantial effect on the 



mea sured extinction law jKenvon. Lada &: Barsoriy 1998). Kenyon. Lada &: Barsony (1998) 
and I Gomez &: KenyonI (120011 ) observed the nearby star-forming regions p Oph and Cha 
I, respectively, and determined the near-infrared reddening law using the same technique, 
but with different photometric systen as, and find a sub stantial difference in the reddening 
law between these two lines of sight. iNaoi et al.l (120061 ) reobserved p Oph and Cha I with 
one photometric system and used the same techniqu e, the same off-cloud regions, an d the 
same limiting magnitud es as the previous studies by iKenyon. Lada fc Barsonyl (119981 ) and 
Gomez fc KenyonI (12001! ) b ut found no variat i on in the near-infrared reddening law between 
these two molecular clouds. iNishiyama et al.l (120061 ) r neasured the extin ction law toward the 
galactic center using the same photometric system as iNaoi et al.l (120061 ) and their results are 
consistent with the reddening law being universal between the Galactic Center and the p 
Oph and Cha I s tar-forming regions . Evidence for a variation in the near-infrared reddening 
law comes from iRacca et al.l (120021) who obser ved the Coalsack Globule 2 using the same 
photometric system as lGomez &: KenyonI (I2OOII ). and used the same technique to determine 
the reddening law, but found a significant difference between the reddening laws along these 
two lines of sight. 

To isolate the infiuence of the near-infrared extinction law on our measured color-excess 
ratios we have measured Eh-KsI Er^-j ^ manner similar to the previous color-excess 
ratios. Since we did not require a J band detection for the Orion A and Ophiuchus regions, 
we excluded these regions. NGC 2024/23, NGC 2068/71 and Serpens have color excess ratios 
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of -0.333 ± 0.004, -0.335 ± 0.003 and -0.328 ± 0.004. The /=284° off-cloud region has a 
color excess ratio of —0.326 ± 0.004. No significant difference exists in the near infrared 
color excess ratios among these lines of sight within uncertainties, indicating that there is no 
significant variation in the near-infrared extinction law through these three molecular clouds 
and through the interstellar medium. 

We examine the mid-infrared portion of the extinction law by forming color-excess 
ratios that include only IRAC bands from the already measured color excess ratios. From 
the definition of the color excess ratio 

R, = = ~ (3) 



Eks-1\] - ^> 



we form the following quantity: 



By choosing Al = 3.6/im, A2 = A3 = 4.5/im and A4 = 5.8/xm we have a measure of the IRAC- 
only color-excess ratio. By combining color-excess ratios rather than using the extinction law 
measurements, Ax/Ak^, we reduce the uncertainty in the IRAC-only color-excess ratio. The 
results are listed in Table H] for all five young clusters as well as for the Z=284° off-cloud line 
of sight. While the uncertainties are large, there is clear evidence for a separation between 
the extinction law measured toward the molecular clouds and the extinction law for only 
the diffuse ISM. The tabulated values of this ratio demonstrate that the variations between 



our color-excess ratios and those of llndebetouw et al.l (120051 ) cannot be explained solely by 



a change in the value of Ah/Ak^, and that the form of the interstellar extinction law must 
differ between the diffuse ISM, Ophiuchus, and the remaining four molecular clouds. It is 
also apparent that the Ophiuchus cloud reddening law differs from the reddening law for the 
other four young clusters. 



4.2. The Extinction Law beyond 8/im 



Extending our analysis out to 24/im we measure Eh-k/ Ek-[24\ = 0.98 ± 0.04 and 
1.14 ± 0.06, A24/AK, = 0.44 ± 0.02 and 0.52 ± 0.03, assuming Ah/Ak, = 1.55, for Serpens 
and NGC 2068/71 respectively. The large uncertainty for both regions' color excess ratios 
reflects the small number of sources used in the fitting procedure. Figure [TD] shows the fit 
to the reddening bands at 24//m, similar to Figures [2] - [61 Figu r e [9] i llustrates that our 
24/xm determinations are in reasonable accord with those of iLutzl ( 119991 ) shown as squares 
in Figure [9l given the large uncertainties, as well as with the measurement at 15/im by 
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Jiang et al.l fl2006h . Ai^/Ak = 0.15 ±0.07, shown as a circle in Figure El iHuard et al.l (120071 ) 
also measured the extinction law at 24/im and found it to be consistent with our results. 



Rosenthal. Bertoldi fc Drapata (120001 ) measured the extinction law from 2 — 30yum to- 



ward OMC-1 using the relative intensities of H2 rotation and ro-vibration lines. They fit a 
functional form assuming a power law for A < 6.5/xm and included the silicate features at 
9.7/im and 18/im, as well as a water ice feature at 3.1/im. The strength of the 9.7/im silicate 
feature was determined from a fit to their observations, while the relative strength of the 
18/im and 9. 7um silicate feature w as taken to be AisM a. j = -44, an average of previou s 



results from iDraine fc Led (119841 ). iPegourie fc Papoularl (119851 ) and IVolk fc Kwokl (119881 ). 
The shape of the 9.7/xm feature was chosen to match the Trapezium emissivity profile, which 
is more appropri ate for extinction d ue to dense molecular clouds rather than the narrower 
/i Cephei profile (IBowey et al.lll998l ). While the assumption of a power law from 3 — 8/ xm is 



clearly not valid, t he extinction law from 10— 30/im, as measured here and by lLuta (119991 ) and 
Jiang et al.l (120061 ). is very well fit by a law do minated by sihcate emission with the adopted 
strengths of the two silicate features given by [Rosenthal. Bertoldi fc Drapata (I2OOOI ) , as can 
be see n in Figure |9] The 9.7Mni silicate feature has been found to vary along different lines of 



sight (IBowey et al. 



19981 ). particularly on the long- wavelength side of the feature, although 



there is insufficient data here to draw any conclusions about the universality of the extinction 
law from 8 — 30/xm. 



5. Conclusion 

We measured the extinction law from 3 — 8/im toward five nearby star-forming regions 
and found it to be relatively internally consistent among different cluster-forming molecular 
clouds although there is evidence that the Ophiuchus extinction law deviates from the other 
four regions at 5.8 and 8/im. The extinction laws we present here d i ffer sy stematically from 
the extinction at IRAC wavelengths derived by llndebetouw et al. J2OO5I ). as seen in both 



the color excess ratios and the relative extinction ratios. To exclude any possibility that 
this difference could be due to the differing methodologies used to derive the extinction law, 
we have also analysed data provided by R. Indebetouw using our methods, and continue 
to find systematic variation. The result could refiect a physical difference in the extinction 
law between the diffuse ISM and molecular clouds, or the differing extinction probes utilized 
in the two studies. Our study of the probes leads us to believe that the major cause of 
the discrepancy is variation in the diffuse ISM vs. molecular cloud extinction laws. The 
observed difference in the extinction law towards Ophiuchus could be caused by diffuse ISM 
behind the molecular cloud that the extinction probes must pass through or differing dust 
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properties between the Ophiuchus cloud and the other four clouds studied. 



The extinction law derived by iLuta (119991 ) agrees more closely with the re-derived 
Indebetouw et al.l (120051 ) extinction law at wavelengths < Sfim . At 3jUm < A < 8/im, 
the opposite is true. At 7/im the relative extinction measured by IJiang et al.l (120061 ) agrees 
with the cluster results. The extension of our extinction law to 24yum was also determined 
from analysis of Spitzer images of Serpens and NGC 2068/71, but further observations at 
this wavelength are required to provide st atisti c ally s ignifi cant conclus i ons. H owever, our 
results at 24/im, alon g with the results of Lutzl (19991) and Jiang et al. (2006), match well 
the functional form of Rosenthal. Bertoldi fc Drapat j (2000) which includes extinction dom- 
inated by silicate features at 9.7 and 18/im. Future measurements of the extinction law for 
10 — 30/im will help to better constrain this fit, as well as to examine any variation among 
different lines of sight. 
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Table 1. Observing Information 



Region 


GTO Program IDs 


Pipeline Version 


Magnitude Zero Points'^ 


Frame Time*^ 


NGC 2024/23 


43 


S9.5.0 


21.93,21.26,19.08,19.44 


12 


NGC 2068/71 


43 


S9.5.0 


21.93,21.26,19.08,19.44 


12 


Orion A 


43 


S9.5.0 


21.93,21.26,19.08,19.44 


12 


Serpens 


6,174,177 


S11.4 


21.99,21.26,19.07,19.44 


12 


Ophiuchus 


6,174,177 


S11.4 


21.99,21.26,19.07,19.44 


12 



For an exposure time of 10.4 seconds and object, inner and outer sky apertures of 2,2 and 



6 pixels respectively. 

'^HDR mode includes the time spent on both the long and short exposure. The effective expo- 
sure time after mosaicing the multiple maps of each region (two for NGC 2024/23, NGC2068/71 
and Orion A, four for Serpens and Ophiuchus) will be longer than this frame time. 



Table 2. Selective Color Excess Ratios 



Re 

NGC 2024/23 1.49 ± 0.02(0.08) 1.17 ± 0.02(0.06) 1.08 ± 0.01(0.06) 1.11 ± 0.02(0.07) 

NGC 2068/71 1.49 ± 0.02(0.07) 1.13 ± 0.01(0.05) 1.03 ± 0.01(0.05) 1.07 ± 0.01(0.05) 

Serpens 1.49 ± 0.02(0.08) 1.17 ± 0.02(0.06) 1.08 ± 0.01(0.05) 1.07 ± 0.01(0.06) 

Orion A 1.51 ± 0.01(0.04) 1.197 ± 0.008(0.03) 1.109 ± 0.007(0.03) 1.116 ± 0.008(0.03) 

Ophiuchus 1.46 ±0.02(0.07) 1.17 ± 0.01(0.05) 1.01 ± 0.01(0.05) 1.01 ± 0.01(0.05) 

Ind05data 1.27 ± 0.02(0.09) 1.11 ± 0.02(0.08) 0.92 ± 0.02(0.06) 0.90 ± 0.01(0.06) 
Ind05paper 1.17 ±0.07 1.0 ±0.03 0.92 ± 0.03 0.92 ±0.04 



Note. — Values "Ind05 paper" are derived from color excess rations from llndebetouw et al.l (120051 ) 
with their uncertainties propogated, and "Ind05 data" are derived from the data R. Indebetouw 
kindly provided. Both are for the / = 284° off-cloud line of sight data. Ind05 data were analyzed in 
the same way that our data were analyzed, in order to enable direct comparison with our results. 
Our uncertainties are for a 68% confidence level with uncertainties for a 99.99% confidence interval 
in parenthesis. 
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Table 3. Relative Extinction 



Region 




^[4.5] /^X, 


^[5.8]/^Xs 


^[8.0] /^X. 


NGC 2024/23 


0.632 ±0.005 


0.53 ±0.01 


0.49 ±0.01 


0.50 ±0.01 


NGC 2068/71 


0.632 ±0.005 


0.51 ±0.01 


0.47 ±0.01 


0.48 ±0.01 


Serpens 


0.630 ±0.005 


0.53 ±0.01 


0.49 ±0.01 


0.49 ±0.01 


Orion A 


0.636 ±0.003 


0.540 ± 0.003 


0.504 ±0.003 


0.506 ±0.003 


Ophiuchus 


0.623 ±0.005 


0.53 ±0.01 


0.45 ±0.01 


0.45 ±0.01 


Ind05 data 


0.57 ±0.01 


0.50 ±0.01 


0.40 ±0.01 


0.39 ±0.01 


Ind05 paper 


0.57 ±0.05 


0.43 ±0.07 


0.41 ±0.07 


0.37 ±0.07 



Note. — Quoted values taken from llndebetouw et al.l (l2005l ) ("Ind05 pa- 
per") are for the / = 284° off-cloud line of sight. The values of A\/Ak^ were 
derived from the measured selective color excess ratios Eh-Ks/Ek^-x ioi the 
five young clusters and the rederived co lor excess ratios for the I = 284° line of 
sight assuming Ah/Ak^ = 1.55 ± 0.08 (llndebetouw et al.ll2005l ). Uncertainties 
on Ax/Ak^ derived here do not include the uii c ertain ty in Ah/A^^- Uncertain- 
ties on Ax/Ak^ taken from llndebetouw et al. J2005h (labeled Ind05 paper) do 
include the uncertainty in Ah/Aks- 
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Table 4. IRAC-only Color Excess Ratio 



Region -E'[3.6]-[4.5]/-E'[4.5]-[5.8] 



NGC 2024/23 


2.58 ±0.66 


NGC 2068/71 


2.49 ±0.56 


Serpens 


2.58 ±0.66 


Orion A 


2.61 ±0.33 


Ophiuchus 


1.25 ±0.13 


Ind05 data 


0.61 ±0.14 



Note. — Color Excess ratio de- 
rived using Eh-kJ Ek,-[\]- See text 
for more details on this calculation. 
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Fig. 1. — J vs. J — Kg diagram for the NGC 2068/71 cluster (plotted in blue), and an 
equivalent area off-field region (black). Red symbols represent cluster stars excluded by the 
H — Kg cutoff. For the off-cloud (black symbols) line of sight, dwarfs occupy the peak at 
J — Ks=0.2, while red clump sources are distributed around J — ^^=0.75. On the other hand, 
toward the clusters (blue) red clump extinction probes are found peaked around 7-/^5=1.05 
at J = 14.5 toward the cluste r. The solid line is a re ddening vector for Ay = 5 extending 
from J = 14.5, J - = 0.75 jRieke fc Lebofskvlll985h . 
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Fig. 2. — 2MASS and IRAC color-color diagrams for NGC 2024/23 non-excess sources. Line 
is the best fit to the data. Sources with H — Kg < 0.2 or Ks — [3.6] < were not included 
in the fit. Slopes are listed in Table [21 
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Fig. 3. — 2MASS and IRAC color-color diagrams for NGC 2068/71 non-excess sources. Line 
is the best fit to the data. Sources with H — Kg < 0.2 or Ks — [3.6] < were not included 
in the fit. Slopes are listed in Table [21 
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Fig. 4. — 2MASS and IRAC color-color diagrams for Serpens non-excess sources. Line is the 
best fit to the data. Sources with H — Kg < 0.2 or Kg — [3.6] < were not included in the 
fit. Slopes are hsted in Table El 
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Fig. 5. — 2MASS and IRAC color-color diagrams for Ophiuchus non-excess sources. Line is 
the best fit to the data. Slopes are listed in Table [2J 
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Fig. 6. — 2MASS and IRAC color-color diagrams for Orion A non-excess sources. Line is 
the best fit to the data. Sources with H — Kg < 0.2 or Kg — [3.6] < were not included in 
the fit. Slopes are listed in Table [2l The dashed line is the boundary used to eliminate the 
branch near H — Kg = 0.25.. 



-34- 



Wavelength (//m) 

10.0 6.67 5.0 4.0 3.33 




0.10 0.15 0.20 0.25 0.30 

1 /Wavelength {fj,m~^) 



Fig. 7. — Selective Color Excess Ratio versus inverse wavelength for NGC 2024/23, NGC 
2068/71, Serpens, Orion A, Opliiuchus and from the / = 284° off-cloud region examined by 
Indebetouw et al.l ((20051), derived using the same method as for the five star- forming regions. 



- 35 - 




Fig. 8. — Residuals {H — Ks)obs — [H — Ks)jnd05 for the non-excess sources in NGC 2068/71 
vs. H — Kg. See text for explanatio n of residuals. The solid horizonal line is flat and is the 
expected slope if the reddening law (jindebetouw et al.ll2005l ) were a perfect fit to the data. 
At larger extinction (larger H — Kg) there is a substantial separation fr om the solid line, 



sugges ting the appropriate reddening law has a different slop e than that of lliidebetouw et al. 



( 120051 ). Vertical dashed line shows the color for a K2 giant (iBessell &: Brettlll988l ). taken to 
be the typical spectral type for a red clump star, extinguished by ^4^=3 using the extinction 
law of iRieke fc Lebofskyi (11985 ). Av=3 is the boundary for ice mantle growth within a 



molecular cloud (IWhittet et al.ll200ll ). 
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Fig. 9.— Extinction in IRAC bands versus wavelength for each of NGC 2023/24, NGC 
2068/71, Serpens, Ophiuchus, Orion A and t he extinction law rederived here for the / = 284° 

(I2005 '). The squares are the extinction law toward 

Lutzl (119991 ) using hydrogen recombination lines. On 
10/im, while on the bottor n is the extinction law 



line of sight from llndebetouw et al 



the galactic center as measured by 
the top is the extinction law from 1 
from 7 — 30/im. The extinction law measured at 7/xm and 15/im by iJiang et al.l (120061 ) are 
labeled with circles in the two plots. On the bottom, the e xtinction law between 7 and 
24/xm is shown for NGC 2068/71, Serpens and as measured by iLutzl (119991 ) with diamonds, 
tri angles and squares respectively. The lin e in the second plot is the extinction curve derived 
by [Rosenthal. Bertoldi fc Drapata (120001 ) using H2 line ratios observed towards OMC-1. 
Horizontal solid lines represent the bandpasses for the IRAC and MIPS filters in the top and 
bottom plots repectively. 
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Fig. 10.— 2MASS and MIPS 24/im color-color diagrams for NGC 2068/71 (top) and Serpens 
(bottom) . Line is the best fit to the data. 



